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Abstract. Limit periodic point sets are aperiodic structures with pure 
^S ■ point diffraction supported on a countably, but not finitely generated Fourier 

T-H ' module that is based on a lattice and certain integer multiples of it. Exam- 

>--^ ■ pies are cut and project sets with p-adic internal spaces. We illustrate this 

^ ' by explicit results for the diffraction measures of two examples with 2-adic 

r- { ' internal spaces. The first and well-known example is the period doubling 

h-^ I sequence in one dimension, which is based on the period doubling substitu- 

lf> • tion rule. The second example is a weighted planar point set that is derived 

from the classic chair tiling in the plane. It can be described as a fixed point 
(-H ' of a block substitution rule. 

-(— > . 

c3 . 1. Introduction 

I '~'i ' The diffraction measure is a characteristic property of a translation bounded 

measure uj on Euchdean space (or on any locally compact Abelian group) and 
'TT ! has important applications in crystallography because it describes the outcome 

r~^ ■ of kinematic diffraction (from X-rays or neutron scattering, say). It is the 

(^ , Fourier transform 7 of the autocorrelation measure 7 of w. Following the dis- 

r^ ■ covery of quasicrystals (which are generally non-periodic but nevertheless show 

^^ ■ pure point or Bragg diffraction), measures with pure point diffraction have been 

^ I investigated thoroughly. In parallel, there has also been progress in understand- 

f— ^ ' ing singular continuous and absolutely continuous diffraction; see [HI H] and ref- 

T-H . erences therein. Let us mention in passing that making a connection with the 

^ I spectra of Schrodinger operators is tempting, but has so far eluded all attempts 

•'~j ■ to substantiate it. 

r^ [ In this article, within the realm of pure point diffractive systems, we consider 

C^ ' two examples of limit periodic structures, the period doubling sequence in one 

dimension (Section [3]) and point sets related to the chair tiling in two dimensions 
(Section [4]). The corresponding point sets arise from fixed points of (block) 
substitution rules, but also have a description as cut and project sets with 2-adic 
internal spaces. We begin with a concise survey of the relevant quantities, where 
we restrict ourselves to measures that are supported on integer lattices. Then, 
the two main examples are presented in an informal way, though all results 
are rigorous (with the underlying mathematical framework being available via 
[UlElll])- A more detailed exposition will appear in [5]. 

2. Autocorrelation and diffraction of Dirac combs 

Consider the (weighted) Dirac comb cj = X^xgZ'* ^(^)'^2^' with 6x the nor- 
malised point (or Dirac) measure at x E Z and (in general, complex) weights 
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w{x). The natural autocorrelation (if it exists) of uj is defined as the Eberlein 
convolution 

[1) ^ = CO ® UJ := hm 



N^ocYol{[-N,NY)' 

which we view as a measure on M'^. Here, w^ = a;|r_^^^id denotes the restriction 
of u to the closed cube [— A^, N^, and Ji is the 'flipped over' measure, defined 
by Jj-ig) = fJ-ig) with g a continuous test function of compact support and 
g{x) = g{—x), the bar denoting complex conjugation. In all cases considered 
below, the limit in ([T|) exists. Due to supp(6t;) C Z,'^, the autocorrelation is then 
of the form 

(2) j = Y,viz)6. 

with the autocorrelation coefficients 

(3) ri(z) = lim — — -j > w(x)w(x — z). 

The existence of the limit in ([1]) is equivalent to the existence of the autocorre- 
lation coefficients 77(2:) for all z G Z . 

By construction, 7 is a positive definite measure, and its Fourier transform 7 
thus exists. The latter is a positive, translation bounded measure on M'^ (by the 
Bochner-Schwartz theorem), which is known as the diffraction measure of cj. It 
corresponds to the (kinematic) diffraction from the measure cj; compare [8j for 
background and |1H [6] for this approach to diffraction theory. The diffraction 
measure possesses the unique decomposition 

(4) 7=(7)pp + (7)3c+(7)3c 

into its pure point, singular continuous and absolutely continuous parts, the 
latter splitting relative to Lebesgue measure, the natural reference measure for 
volume in Euclidean space. Lattice periodic measures in R'^ give rise to pure 
point diffraction, as do regular model sets with Euclidean internal space M"^, 
where the diffraction measure is supported on a Z-module of finite rank ?n, + d; 
see [6] and references therein for more. The examples discussed below are also 
based on model sets, but with 2-adic internal spaces; they are limit •periodic^ 
and thus exhibit pure point diffraction where the corresponding Fourier module 
is countably, but not finitely generated. For both examples, the explicit form 
of the diffraction can be obtained either from the cut and project description 
[6] or from the inflation structure [lOj; in this note, we use the latter approach. 

3. The period doubling sequence 
The "period doubling sequence is defined by the primitive substitution rule 



(5) g: 



b t-^ aa 



on the two- letter alphabet {a, b}. It is related to the classic Thue-Morse system 
(as a factor under a 2 : 1 mapping; compare [3] and references therein) and 
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defines a strictly ergodic dynamical system via the orbit closure of a suitable 
bi-infinite sequence. The latter can be obtained from a fixed point of q^ via the 
iteration 

(6) a\a — > abaa\abaa — > . . . — > w = g (w), 

with convergence in the (obvious) product topology. Here, | marks the origin 
of a two-sided sequence w = . . . w_2W_i\wqWi . . ., and we start with the legal 
seed a\a. The half-infinite word v = WqWi ... is the (unique) one-sided fixed 
point of Q. 

To the sequence w, we associate two point sets yl^ C Z with i € {a, b} via 

Ae = {neZ\wn = i}, 

so AaUAf) = Z, where U denotes the disjoint union of sets. The geometric fixed 
point equation for g^ now implies the equations 

(7a) Aa = 4Aa U {4Aa + 2)0 {4Aa + 3) U AAb U {4Ab + 2) , 

(7b) ylfe = (4yl„ + l)U(4A + l)U(4ylfe + 3), 

which can be decoupled and solved by iteration, yielding 

(8) A = {-l}u|J^^^(2.4^Z + (4*-l)), A = |J^^^(4^Z+ (2.4^-1-1)), 

where the singleton set {—1} has to be added to the fixed point (to either Aa 
or A),), as it is only a limit point under 2-adic completion [7J. 

We attach a Dirac comb to Ag by u^ := 5ai = YlneA ^n ^^ ^ ^ W^b}. A 
general weighted Dirac comb with weights a (for letter a) and /3 (for letter b) 
can then be expressed as a; = aw^j + (3u)f^. Its autocorrelation measure exists 
and has the form of Eq. ^ (with d = 1) with autocorrelation coefficients ri{z) 
according to ([3]), where the weight w{x) € {a,/3} is chosen according to the 
letter at position x € Z. 

For simplicity, we work with the one-sided fixed point v (which has the same 
autocorrelation). We first consider weights a = 1 and /? = — 1, and denote the 
autocorrelation coefficients for this balanced case by rj^{z). The correspond- 
ing sequence {fjjjgpj of weights Vi € {±1} satisfies the recursions V2n = 1 
and 1^2^+1 ~ ~''^n f°^ n > 0. This implies a recursion for the autocorrelation 
coefficients, 

(9) r/_|_(2r?T,) = -(l + r7_|_(r7T,)) and ri_^{2m + 1) = , 

for all m > 0. In particular, 'i]^(0) = 1, and we also have i]^{—m) = r/_|_(m). 
When m = 2^' (2s + 1) with r, s > 0, the recursion leads to 

1 

ri_^{m) = 1 



3 • 2'^-2 ' 

which is independent of s. Going back to general weights a and f3 amounts to 
considering the bi-infinite sequence ^((a + /3)1 + (a — (3)w), or alternatively 
^((a + /3)1 + [a — /3)u) for its one-sided counterpart, which results in 

(10) r? = i(Q-/3)V + n(3(« + /3)' + 2(a'-/32))l. 
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Figure 1. Sketch of the diffraction amphtudes |^(A;)| from 
Eq. (J12p for the period doubhng chain, with k € [0,1]. The 
diffraction pattern repeats Z-periodically along the real line [2]. 



Taking the Fourier transform, the diffraction measure is obtained as 

|2 



(11) 



7= Y.\aAik)+/3B{k)\U,„ 

k&L® 



where L® = [j^^^Z/2' = {p | (r = 0, m G 
Fourier module of the period doubling sequence 

2 



Z) or (r > 1, m odd)} is the 
The amplitudes read 



(12) 



A{k) 



Jlirik 



-2y 



and B{k) = 5^Q- A{k) 



where we implicitly refer to the parametrisation of L® . 

To derive these formulas, one observes that a finite union in Eq. ([8|) would still 
be a periodic set, whose diffraction can be calculated by direct Fourier transform 
(which is then a periodic point measure) followed by taking the absolute square 
of the amplitudes (note that the Wiener diagram is still commutative in this 
case). This results in a sequence of pure point measures that converge not 
only in the vague topology, but also in the stronger norm topology towards the 
diffraction measure in pT|) . with the amplitudes as in p^ : see |6| for details. 
Since pure point measures are closed in the norm topology, this also gives the 
constructive proof that L® is indeed the Fourier module. 

As mentioned above, the same result can also be derived via the description 
of the periodic doubling sequence as a model set with 2-adic internal space; see 
[71 E] for a proof. A sketch of the diffraction pattern is shown in Figure [TJ 

The presence of the powers of 2~* in the Fourier module shows the limit 
periodic structure; since s can be arbitrarily large, the Fourier module cannot 
be finitely generated. Nevertheless, the diffraction measure is pure point, with 
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amplitudes that reflect the inflation symmetry with integer multiplier 2, as 
visible in Figure [H 



4. Block substitution for the chair tiling 
Let us move on to our planar example. Consider the block substitution rule 



(13) 



^ 



1 
3 



1 ^ 



1 2 
1 



2 ^ 



1 2 

2 3 



3 ^ 



3 2 
3 



on the alphabet {0, 1, 2, 3}. Iterating the rule, starting from a legal seed, 



(14) 



3 2 
3 


1 
3 


1 2 

2 3 


1 2 
1 



3 2 12 


12 10 


3 2 3 


10 3 


10 3 2 


10 3 2 


3 3 


3 3 


12 12 


12 12 


12 3 


12 3 


12 3 2 


10 12 


2 3 3 


3 1 



produces a sequence of growing square-shaped blocks (with lines denoting the 
coordinate axes, so the origin is located in the centre of the blocks), where 
each block reappears in the centre of the next. This sequence thus converges 
in the product topology, where the limit is a fixed point of the block substitu- 
tion rule that covers the entire square lattice. This fixed point has a D4 colour 
symmetry: An anti-clockwise rotation through tt/2 corresponds to the cyclic 
permutation (0123), while a reflection in the central horizontal gives the per- 
mutation (01) (23). The pattern is thus invariant under any L)4 transformation 
followed by an appropriate permutation. 

If we choose a representation with four different unit squares and take the 
lower left corners as their reference points, the fixed point defined by (J14p results 
in a partition Z^ = ylQUyl^^LJylgUylg of Z-^ into four point sets of equal density 
i. As a result of (flBj) . they satisfy the fixed point equations 

(15a) Aq = 2Aq\J {2AQ + u + v)U2A^\J2Ar^, 

(15b) A^ = {2Aq + v)0 {2A^ + n) U {2A^ +v)0 {2A^ + v), 

(15c) A2 = {2A-^ + n + v) U 2/12 U (2^2 + u + v)\J {2A^ +u + v), 

(15d) ylg = {2Aq + n) U (2yl2 + u)(j {2A.^ + n) U {2A^ + v), 



^ and V 

T 



(0,1: 



with u = (1,0) 

I? as r+ = {(a:i,X2J^ ^ £," \ xi 

1? = r+ur_ , 

decoupled equations 



T 



as well as ylQUyl2 



Denoting the even and odd sublattices of 
X2 = mod 2} and /^_ = n + I^, we have 
= Fj^ and yl^Uyl3 = r_. This leads to the 



(16) Aq = 2r_ U (2ylo + {0, n + v}) and A^ = (2r+ + v) U {2A^ + {u, v}) , 
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together with A2 = i~Y \ ^0 ^^^ ^3 ~ ^- \ ^i- ^^ exphcit solution is (with 
Sx := {sx\se S} for sets S C Z and Sr := {0, 1, . . . , 2'"- 1} for r > 0) 

(17a) Aq = noiu + v)[j\J^^^{2'+^r_ + Sriu + v)), 

(17b) Ai+v= No(u-?;)u|J {2''-^^r_ + Sr{u - v)), 

(17c) A^ + u + v =-Nq{u + v)0\J {2''+^r_-Sr{u + v)), 

(17d) As + u= -No(n - «) U |J,,>o (2"^'^- " 5r(n - v)), 

where we used /+ + u = /+ + f = /^_ = F^ + u + v. Each of the infinite 
unions on the right-hand side defines a point set of density j. The sets of 
points (of density 0) along the diagonals have to be added for our fixed point 
(|14p . because they are not contained in any of the lattice translates (but are 
covered by the 2-adic completion of the infinite unions, similar to the situation 
of the period doubling sequence). The solutions obey the relations yloUyl2 = F^ 
and ^iLJylg = r_. 

Now, we are going to calculate the diffraction measure of the Dirac comb 

3 

(18) U = Y,»^h. 

j=0 

with the four point sets A- from Eq. ()17p . and arbitrary complex numbers aj. 
Define Qr{x) := 2''+ir_ + SrX for r € Nq and x G Z^. Using F^ = ^F^ 
(where the superscript * denotes the dual lattice), we calculate 

2''-l _ 2^+iTrikx ^ ^ 

? _ V^ -2nimkx Is: \^ _ ^ "^ -2^+'^mku r 

°Qr{x) ~ Z^ ® \P2^+^(r++u)) - 'Y_ g-27rifex ® °2^+^r+ 

m=0 
]^ _ ^-2^+^Tv\kx g-2''+27rifcM 

(1^) ^ 1 _ Q-2TTikx ^27+3 '^r+/2'-+2' 

where we again write density factors as functions of k. The first fraction eval- 
uates as 2^ whenever fcx € Z (by an application of the I'Hospital rule). The 
resulting Fourier module is 

(20) 

^®= U^ = ^'^U,>i{^("''^) |m,nGZwith gcd(m,n, 2) = l}, 

s>2 

which is the coarsest group that contains all positions where the Fourier trans- 
form in (|19p has a point measure. 

The diffraction measure of our weighted Dirac comb uj is of the form 

3 

(21) € = E|E«*^*(^) 



fcgL® i=0 



2 
Sk, 



where Ai{k) is the amplitude that belongs to Ai. Its derivation is analogous 
to the method explained above for the period doubling sequence (via periodic 
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approximants that converge in the norm topology). This method also provides 
a proof for the correctness of (|20l) . 

The amplitudes themselves can be calculated in the same way as for the 
period doubling sequence above. For any k € L®, there is a minimal s > such 
that k E r'+/2*+-^, and then k € /+/2^'+^ for all r > s. The amplitudes result 
from summing the corresponding contributions from Eq. (J19p . with x = u -\- v 
for i = 0, x = u — V for i = 1, x = —u — v for i = 2 and x = —u + v for 
i = 3. In addition, the results have to be multiplied by an overall phase factor 
that reflects the translation vectors on the left-hand sides of the four relations 
in Eq. ()17p . We omit the detailed (and somewhat lengthy) calculations and 
simply state the result. 

To this end, the integer and half-integer points of L® are better treated 
separately. Whenever A: G Z^, 



(22) 



Ao{k) = Ai{k) = A2{k) = A^ik) = ^. 



If /c € 2-f+ \ ^^) one has k = i^irn, n) with m and n odd, so that kx G Z for all 
X G r^. This gives 



Ao{k) = A2{k) = - and Ai{k) = As{k) 



(23) 

while k = 2 (?^) n) with m + n odd results in 



(24) Aoik) 



A2ik) 



Aiik) 



(-ir 



and ^^{k) 



{-ly 



In particular, ylo(^) + ^2(^) = 0, and also Ai{k) -|-^3(A;) = 0, as m -|- n is odd. 
The remaining elements of L® are of the form k = (m, n)/2^ with s > 2 and 
m,n £Z subject to the condition gcd(m,n, 2) = 1. The amplitudes A^i^k) and 
A2{k) depend on m + n. They satisfy A2{k) = —AQ(k) and 



(25) Ao{k) 



0, 



1 



2 1 

'4^~ 



(_l)(m+n)/2 



1 



m+n 



li" 1 



£. 



m+n ' 



if 2^* I (m + n), 

if m + n even with 2^ f (m + n) , 

if m + n odd. 



where e^ = exp(— 27ri/2*) is a root of unity. Similarly, one obtains the relation 
Asik) = —Ai{k) together with 

'o, 



(26) Ai{k) 



In particular, Oq 



1 



2 1 

4^~ 



(_l)(m-")/2 



1 



m—n 



U'^ 1 



m—n ■ 



if 2* j (m — n), 

if m — n even with 2* f (m — n) 

if m — n odd. 



Qi 



an 



03 = 1 leads to 7^ = 5^2 , as required. Moreover, 



setting Oq = 02 = 1 and a^^ = Q3 = gives the diffraction 7^^ = 5 



while the alternative choice a 







Or, 



together with a^^ 



r+ 



Oq 



(using the parametrisation of k as above) 7^ 



-2TTik 



■T- 



™5 



^+" 



(-1)- 



2"r+/2' 
1 yields 

V+/2- 
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Figure 2. Illustration of the diffraction measure for the Dirac 
comb ([2T|) with weights Oj = P , where j G {0, 1, 2, 3}. A peak at 
k £ [—1, 1]^ is represented by a black disc, centred at k, with an 
area proportional to the peak intensity. The complete pattern 
is i/V-periodic. 



An interesting case results from the balanced choice a ■ = i^ , which produces 
extinctions for all k E 2-^+' '^^^ corresponding diffraction measure is illustrated 
in Figure [21 It is Z)4-symmetric, which is not the case in general. 

For a generic choice of the weights, our Dirac comb (jlSp is mutually locally 
derivable (MLD) with the chair tiling. As it is based on four subsets of I?, the 
diffraction measure 'y^ is always Z^-periodic [2]. This is the generic maximal 
translation symmetry. However, since yloUyl2 = /+, the diffraction of any 
weighted Dirac comb with only these two components is periodic under the 
dual lattice F^ = ^r_^_. Likewise, the corresponding statement holds for the 
diffraction of weighted Dirac combs with components Ai and A3. This is also 
the translation symmetry in Figure [21 which is a consequence of the particular 
choice of the weights here. 

Any choice of the weights refers to a discrete structure that can be viewed as 
a decoration of the chair tiling, and is thus locally derivable from it. Conversely, 
a typical decoration of the chair tiling will be locally derivable from our coloured 
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point set, but need not be realisable as a simple decoration of the square-shaped 
building blocks. In this sense, our above results constitute only a first step of 
the complete diffraction analysis of the chair tiling. 

5. Summary and Outlook 

Limit periodic structures form an interesting class of pure point diffractive 
systems. Explicit results for the diffraction can be obtained, at least for ex- 
amples such as those discussed in this paper; related examples can be found in 

mm- 

Limit periodic tilings also arise in the context of Wang tiles. Raphael Robin- 
son's tilings consisting of six square-shaped prototiles [H], up to rotations and 
reflections, is one example. In this case, the limit periodic structure is apparent 
from a hierarchical pattern of interlocking squares of larger and larger sizes. 
Recently, a similar tiling with a single hexagonal prototile, up to rotations and 
reflections, was discovered |151 I16j . with a hierarchical pattern of interlocking 
triangles. It is related to the half-hex inflation that was studied in detail in ^, 
where also its limit-periodic structure was derived, and to Penrose's (l-|-e + e^) 
tiling [12]. 

The new example is particularly intriguing because the hexagon is the first 
example of a single aperiodic tile (provided you count the tile and its mirror 
image as a single tile), in the sense that there exist local rules that allow the 
existence of tilings covering the entire plane, but no such tiling can have any 
non-trivial period. The local rules involve not just neighbours, but also the 
next corona, and thus cannot be encapsulated in the shape of a tile whose 
interior is connected; see |15] for details. It would be interesting to study 
the diffraction properties of these tilings, which will refiect the limit periodic 
structure apparent in the patterns of squares and triangles. 
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